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Abstract 


The Advanced Very High Resolution (AVHRR) Polar Pathfinder Data (APP) provides 
the first long time series of consistent, calibrated surface albedo and surface temperature 
data for the polar regions. Validations of these products have consisted of individual 
studies that analyzed algorithm performance for limited regions and or time periods. This 
paper reports on comparisons made between the APP-derived surface albedo and that 
measured at fourteen automatic weather stations (AWS) around the Greenland ice sheet 
from January 1997 to August 1998. Results show that satellite-derived surface albedo 
values are on average 10% less than those measured by the AWS stations. However, the 
station measurements tend to be biased high by about 4% and thus the differences in 
absolute albedo may be less (e g. 6%). In regions of the ice sheet where the albedo 
variability is small, such as the dry snow facies, the APP albedo uncertainty exceeds the 
natural variability. Further work is needed to improve the absolute accuracy of the APP- 
derived surface albedo. Even so, the data provide temporally and spatially consistent 
estimates of the Greenland ice sheet albedo. 


? 


(1) Introduction 


Surface albedo, the hemispheric reflectivity integrated over the entire solar 
spectrum, is a fundamental component for determining the Earth’s climate. It is a 
parameter needed in both global and regional climate models, and for determining 
surface energy balance. In the polar regions, the high surface albedo acts to cool the 
lower troposphere and therefore acts as an important boundary condition for energy and 
heat exchange between the mid-latitudes and the polar regions. The high surface albedo 
of snow and ice-covered surfaces allows for very little solar energy to be absorbed by the 
snowpack. However, as the snowpack begins to melt the reduction in surface albedo 
associated with that melt allows significantly more solar radiation to be absorbed. 
Consequently, snowmelt comprises an unstable positive feedback component of the 
climate system, which amplifies small perturbations to that system. 

Current estimates of snow albedo fall short of the accuracy needed for climate 
studies. The accuracy of snow albedo in General Circulation Models (GCMs) is still 
found to be poor because of inadequate parameterizations ( Roesch et al., 1999). Since 
climatic warming is closely connected with alterations in the surface energy balance and 
hence, surface albedo, it is essential to monitor these changes in response to changing 


climatological conditions as accurately as possible. 


Although the surface albedo has been routinely observed for some time, maps of 


surface albedo in the polar regions have not been routinely provided. One of the AVHRR 
Polar Pathfinder (APP) products are gridded maps of the surface albedo for the polar 
regions spanning 1981 to 1998 based on observations from the Advanced Very High 
Resolution Radiometer (AVHRR) on board NOAA POES satellites. This NOAA/NASA 
Pathfinder program project is currently generating multi-year time series of moderate- 
resolution products derived primarily from AVHRR data {Meier et al., 1996; Hutchinson 
andScambos, 1997 ; Maslanik et al., 1997). 

One of the Polar Pathfinder projects' main goals is to provide consistent data sets 
for the polar regions. The Polar Pathfinder Program helps ensure consistency between 
current Polar Pathfinder products and provides a transition to similar products derived 
from existing and forthcoming satellites, such as from the Moderate-Resolution Imaging 
Spectroradiometer (MODIS). To meet this goal requires validation of the existing Polar 
Pathfinder Products using field observations. 

Validations to date for the APP surface albedo consist of individual studies that 
analyzed algorithm performance applied to AVHRR data for limited regions and/or time 
periods (e g. Meier et al., 1996; Maslanik et al., 1997). While the surface albedo maps 
over the ice sheets appear qualitatively correct, they are of poorer absolute accuracy than 
was found in previous studies that examined albedos of sea ice and snow on land. 
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This paper reports on comparisons made between the 1.25 km AVHRR Polar 


Pathfinder-derived surface albedo and that measured at fourteen automatic weather 
stations on the Greenland ice sheet (Figure 1) from 1 January 1997 to 6 August 1998. 
The main purpose of the study is to validate the existing APP-derived surface albedo 
estimates, but results of the comparison will also serve to improve AVHRR-derived 
surface albedo maps for the Greenland ice sheet. 

(2) Data 

2.1 Satellite Data 

The AVHRR Polar Pathfinder data set includes multichannel data and 
geophysical products derived from AVHRR flown on the NOAA/TIROS operational 
meteorological satellites (NOAA-7, -9, -1 1, and -14). The APP product provides 
continuous daily data from July 1981 through summer 1998, with twice per day surface 
temperature and albedo and daily sea ice velocities for both polar regions. Considerable 
effort has been made to provide a consistent data record by cross-calibrating individual 
instruments and minimizing processing changes over the July 1981 to August 1998 
record. The APP products were processed at the Colorado Center for Astrodynamics 
Research (CCAR) at the University of Colorado and are being archived at the National 


Snow and Ice Data Center (NSIDC). 



The APP data set used in this comparison study is the Version5 1 ,25-km Polar 
Pathfinder Product data, derived from the 1-km NOAA-14 AVHRR Level lb scenes. 
The Version5 product for the Northern Hemisphere contains 2 Granules per day, near 
0400 and 1400 Local Solar Time (LST). Each Granule contains the following products: 
broadband surface albedo, surface skin temperature, AVHRR channels 1-5, solar zenith 
angle, relative azimuth angle, satellite elevation angle, cloud mask, orbit mask, ancillary 
information file and universal time coordinate. The data are gridded onto a Lambert 
Equal-Area 1.25-km grid derived from SSM/I Pathfinder's EASE Grid {Armstrong and 
Brodzik, 1995). Northern Hemisphere grid coverage extends from the North Pole to 
50°N. Data at 1.25-km resolution are currently available from 1993 to 1998. 

AVHRR channels 1 and 2 are calibrated using the most recent coefficients 
developed by the NOAA/NASA AVHRR Pathfinder Calibration Working Group {Rao 
and Chen, 1999). The data are then georeferenced to earth coordinates using an orbital 
ephemeris model with orbit and clock time corrections {Baldwin and Emery, 1995). 

APP surface albedo is derived from top-of-the-atmosphere (TOA) AVHRR 
channel 1 and 2 radiances using several steps as depicted in Figure 2. First, the channel 1 
and 2 calibrated TOA reflectances are converted to a broadband reflectance based on a 
linear relationship between the channel 1 and 2 reflectances and a TOA broadband albedo 
for snow-covered surfaces. To develop the regression equation, the STREAMER 
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radiative transfer model (Key and Schweiger, 1998) was used to simulate the TO A 
reflectances over a broad range of viewing and illumination angles and atmospheric 
conditions (aerosol optical depth and water vapor). 

Next, the directionally dependent TOA broadband albedo is converted to a 
hemispheric albedo using ERBE/GOES-derived anisotropic correction factors (Suttles et 
al., 1 988). The work by Suttles et al. ( 1 988) updates earlier work by Taylor and Stowe 
(1984) to convert the directional reflectance to hemispheric albedo. A tri-linear 
interpolation routine was used to derive the conversion factors as a function of satellite 
zenith angle, solar zenith angle and relative azimuth angle. Former to this technique, the 
Lindsay and Rothrock (1994) parameterization of the Taylor and Stowe (1984) 
anisotropic correction factors for snow surfaces was used. The two methods are 
compared later. 


Finally, the surface broadband albedo is calculated by correcting for atmospheric 
extinction after Koepke (1989). 


Surface albedo = (TOA broadband albedo - a) / b 


Eq. 1 


The atmospheric correction coefficients 'a' and 'b' were developed from regression of 
output from the STREAMER radiative transfer model using a variety of atmospheric 
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conditions and viewing geometries. The atmospheric correction is a function of solar 


zenith angle, water vapor and aerosol optical depth. 

In the atmospheric correction, the aerosol optical depth and ozone concentration 
are prescribed from climatological averages to be 0.06 and 0.24 cm, respectively. Since 
no atmospheric profiles of water vapor are available for the Greenland ice sheet, 
estimates of cleaf-sky total precipitable water vapor are computed based on differences in 
temperatures between channels 4 and 5. The conversion from brightness temperature 
difference to precipitable water vapor was derived from regression of Lowtran-modelled 
brightness temperatures at 1 1 and 12 jam and precipitable water obtained from Arctic 
rawinsonde data. No rawinsonde data from Greenland were used in deriving the 
conversion. 

It is important to note that the APP surface albedo is an “apparent” surface 
albedo, the ratio of upwelling irradiance to downward irradiance at a specific solar zenith 
angle and wavelength range. It is the same albedo that is measured by radiometers at the 
surface and is a function of atmospheric conditions and solar zenith angle. 

Cloud screening is applied to the data using algorithms from the Cloud and 
Surface Parameter (CASPR) toolkit (Key, 1999) However, note that the same surface 
albedo algorithm is performed on all pixels, regardless of cloud cover. Instead of 
computing a “cloudy” sky albedo, an accompanying cloud mask based on CASPR is 
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2.2 Station Data 

Observations of incomi ng and reflected shortwave radiation (0.4-1. 1 pm) are 
available from The Greenland flimate Network (GC-Net) ( Steffen ( * al ., 199 The G - 
Net is part of the Program for Arctic Regional Climate Assessment (PARCA) (Thomas et 
al. 2000 submitted to JGR May 3 1 2000) to provide ground truth information for existing 
and forthcoming satellites in addition to surface energy and mass balance measurements 

There are currently 18 operational GC-Net sites, including the 3 stations installed 
during the spring of 1999 (see T die 1 and Figure 1). The stations contain instruments to 
measure the surface radiation balance, turbulent fluxes, conductive heat flux in the snow 
pack, profiles of temperature, humidity and wind speed and direction, snow height and 
pressure. Hourly averages are tra smitted via a satellite link (GOES c, ARGOS) 
throughout the year. In addition, e measurements are stored locally ; solid s 


memory 



sources of variance. However, the overall difference in the average co-registered albedo 


was 0.04. 

Other errors in the measured surface albedo may result from the following: (1) 
imperfect or changing instrument level, (2) the instrument ‘cosine response’ and (3) other 
problems associated with instrument obstruction by frost or reflections/shadows onto the 
instrument from the tower The latter sources of error are infrequent since the 
instruments are mounted at the end of a boom arm pointed south to minimize shadowing 
near midnight However, some problems are still evident in the data, such as rime frost 
on up-looking pyranometers at the NASA-E site from days 104-117 and at the North 
Greenland Ice Project (NGRIP) site from days 61-91 during 1998. 

The primary source of error in the measurement of surface albedo is instrument 
level. The bias in the resultant surface albedo may be amplified if either or both LI- 
CORs are not perfectly level. Such a bias is very difficult to correct. Instrument leveling 
may drift because of snow compaction effects on tower anchoring, and settling of the 
boom arm through time, despite the fact that the boom arm has a cable mounted near the 
tower to provide counter tension. Leveling errors are corrected at each visit (every 1-2 
years). 

Uncertainties in LI-COR measurements increase at low signal levels, which 
typically occur at high solar zenith angles. The uncertainty occurs because the 
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pyranometer response is normalized to a single solar zenith angle, usually 50°. 
Pyranometer response is less accurate for large solar zenith angles. Departures from 
highly accurate instruments in excess of 5% are observed for solar zenith angles greater 
than 75 °. For this reason, we limit our quantitative comparisons to solar zenith angles 
less than 75°. 


(3) Methods 

Because the current method for retrieval of surface albedo is valid for clear-sky 
cases only, accurate cloud detection is critical. The success of cloud detection over snow 
remains a large source of uncertainty in the current estimates of surface albedo using 
AVHRR data. Because uncertainties remain in the APP cloud detection method over 
Greenland, the cloud detection scheme of Box (1997) was used instead. This method is 
based on comparisons between observed incoming solar irradiance with clear-sky 
irradiance calculations from the FluxNet radiative transfer model (Key and Schweiger , 
1998), providing an index of effective cloud opacity. The “effective cloud transmission” 
(t £ ) is defined as the ratio of the measured incoming solar radiation to that computed by 
the radiative transfer model. 


Eq. 2 


T=SW, n . obs /SW 


m.clr 
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A T e value of 1 would imply clear-sky conditions. However, diffuse sky 
irradiance will act to underestimate true cloud amount. Nonetheless, much effective 
cloudiness will be captured since clouds often reduce incoming solar radiation. 

Comparisons with the PSP have indicated FluxNet computed incoming solar 
radiation values greater than the measurements by up to 7%. Given the uncertainty of 
(X) in T efr , T efl >0.8 is used as a binary indication of clear skies. 

(4) Results 

In Figure 3 (a)-(d) the surface albedo at four different AWS sites during the 
afternoon pass in 1997 are shown. The 4 stations, in order of increasing elevation, are (1) 
JAR (962m), (2) Swiss Camp (1 149m), (3) Crawford Point-CPl (1958m) and (4) 
Summit/GISP2 (3 1 50m). These 4 stations represent a transect across the mid-western 
part of the ice sheet (see Figure 1). The Summit and CPI stations are in the accumulation 
facies of the ice sheet. Net annual accumulation at the two sites is 62 cm and 134 cm, 
respectively. The Swiss Camp is located at the mean equilibrium line altitude of the 
Greenland ice sheet and thus experiences a zero net annual accumulation. The JAR site 
is located in the ablation region, and loses approximately 104 cm of snow/ice per year. 
During 1997, the JAR AWS site was situated in a lake of approximately 500m diameter 
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that formed during summer. In 1997, the melt pond filled by 67 cm in 94 hours, 
decreasing the surface albedo from 0.8 when the glacial ice was snow covered, to 0. 14 at 
the time of maximum melt pond depth of —80 cm in 1997. The AWS tower has since 
moved upslope of the lake. The JAR AWS is the only station with a net positive energy 
balance. 

Representative surface albedos are shown in Figure 3 for all coincident 
measurements where both the satellite and in situ data were available and where the solar 
zenith angle was less than 75°. In the Figure we plot the albedo as given by the daily 
surface albedo data files. Thus, no prior cloud filtering was performed in Figure 3. 
Although quantitative comparisons can only be made for “clear-sky” cases, these Figures 
provide a qualitative assessment of how well the APP data represent the general seasonal 
trends in the surface albedo at the different stations. 

The data show that the APP-derived surface albedos produce the same seasonal 
trends in surface albedo as seen in the station measurements. For example, at the two 
sites where substantial summer melting occur (e g. JAR and Swiss Camp) the decrease in 
albedo associated with that melting is also detected by the satellite surface albedo 
estimates. At the JAR site, the APP-derived surface albedo tracks the in situ albedo well, 
showing a steady decline in surface albedo as the surface melts and the melt lake fills up 
with water. On day 176 there is a large drop in the measured surface albedo as the lake 
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reaches its maximum depth. This albedo reduction is detected by the APP-derived 
albedo, albeit not as great as seen by the station. One likely reason that the albedo 
reduction is not as large in the satellite measurement is because of the large satellite 
footprint size. Thus, the satellite measurement contains information about the 
surrounding ice in addition to the surface lake. Albedo increases associated with new 
snowfall events are also detected with the AVHRR measurement, such as on days 215 
and 238. 

Similarly, at Swiss Camp, albedo variations associated with melt and snow events 
are detected in the AVHRR measurements. However, during the beginning part of the 
year in 1997 (March and April) large differences between AVHRR and in situ albedos 
are found. The station measurements show lower than expected surface albedos during 
March and April, (0.70 to 0.75). The GC-Net measurements suggest a surface that has 
undergone melting during winter. Occasionally, mid-winter air temperatures exceed the 
melting point associated with warm air advection in snow and even rain events. There 
was enough snow accumulation since the previous melt season at this station that the 
surface should have been indicative of dry snow. However, analysis of the station data 
suggest that the new snow deposited during March and early April was subsequently 
blown away, leaving a darker, crusted surface exposed. The surface albedos return to 
more “normal” values in April 1997 after a 20-cm snowfall that remained at the surface 
The reason the satellite detects overall greater surface albedos may be a result of the 
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differences in footprint sizes between the satellite observation and the ground 
measurement. Therefore, it’s possible that the satellite might be viewing not only places 
where the new snow has been blown off, but also areas where the snow accumulated 

At higher elevations (~2000m and greater), little summer melting occurs and the 
in situ measured surface albedo does not vary by more than a few percent over the entire 
season. The APP-albedo exhibits similar stability, except for what appears to be about a 
7 to 10-day oscillation in the albedo. This fluctuation in albedo also occurs at the lower 
elevations, but is more readily detectable at the higher elevation stations where albedos 
are less variable. The oscillation is believed to be associated with the variation in the 
satellite elevation angle. For example, in Figure 4 we show the APP-derived surface 
albedo at the Summit together with the satellite elevation angle at this site. The time 
period from approximately day 1 10 to 220 was relatively cloud free, and the in situ 
measurements suggest a relatively constant surface albedo near 0.90. The APP data on 
the other hand show fluctuations in albedo on the order of 7 to 12% (absolute) that appear 
to be associated with changes in the satellite viewing angle. Specifically, a large increase 
in albedo is found when the satellite elevation angle jumps from near 40° (50° off-nadir 
viewing) to 80-90° (near nadir viewing). This suggests that some angular dependency 
remain in the APP-derived surface albedo. 
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In theory, if the conversion is done correctly, there should be no dependence of 
the hemispheric albedo on satellite viewing angle or relative azimuth angle, provided the 
scene is stable. One simple way to test the validity of the angular model is to see if the 
derived albedo changes with viewing angle and relative azimuth angle. In Figure 5 the 
dependence of the APP-derived surface albedo on relative azimuth angle and satellite 
elevation angle is shown using two different angular models: (1) Lindsay and Rothrock 
(1994) and (2) Suttles et al. (1988). In constructing Figure 5, all clear-sky data from the 
Summit, NGREP and NASA-U during 1997 were used. These 3 stations were chosen 
since they are all relatively close to each other and are at high elevations, and thus 
dependencies that might arise from changing surface conditions (e g. melt) are 
minimized. The different symbols in Figure 5 correspond to 3 different relative azimuth 
angle bins: (*) 0-60°, (+) 60-120°, and (A) 120-180°. Discrepancies between the symbols 
reveal to what extent the angular model is not able to represent the albedo dependence 
with relative azimuth angle. Results are shown for only one solar zenith angle bin, 60- 
70°. 


To gain a more quantitative measure of the extent to which the two albedos are 
related we list the clear-sky correlation coefficients computed between the in situ and 
APP-derived surface albedo for all stations during both 1997 and 1998, and for both 
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passes in Table 2. Results are not presented for stations 4, 10, 1 1 and 15, where too few 
overlaps with the satellite data were available and/or there were in situ data problems, 
such as frost on sensors and power outage. Correlation coefficients are listed for clear-sky 
cases only, cases where T efr > 0.8 and where the solar zenith angle is less than 75°. The 
number of observations used in the regression analysis is given in parenthesis. 

In general, the correlation between the two albedos is weak, and at some GC-Net 
sites no significant correlation exists with the APP albedos. This is especially true during 
the morning pass, reaffirming that problems in the albedo retrieval are more pronounced 
at high solar zenith angles. However, this may not be the case since there are less “good” 
data to work with during the early morning pass because we are limiting the analysis to 
solar zenith angles less than 75°. 

Good agreement, with statistical significance at the 0.0005 level is found at sites 
where distinct seasonal variations in albedo occur. These are the sites where summer 
melting occurs; such as at the JAR and Swiss Camp sites. For example, the correlation 
coefficient is 0.87 at the Swiss Camp in 1998. At this station significant melt occurred in 
1998, with 80 cm of snow and 1.4 m of ice melted, resulting in large drops in the surface 
albedo. Less melt occurred during 1997, and a much lower correlation (0.55) is found. 
However, the overall lower correlation coefficient is partially a result of differences 
observed in the early part of the year. Better correlation is found from May through 
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September, but still less than observed during 1998. Statistical significance (based on 
students t-test) is also found at Crawford Point-CPl during the afternoon pass for both 
years and the morning pass during 1997, and at Crawford Point-CP2 and DYE -2 for the 
afternoon pass during 1998. These sites also experienced some slight melting during 
1997 and 1998. Statistically significant correlation was also found at the TUNU-N site 
for the afternoon passes in 1997, a site with little change in surface albedo. 

Figure 6 illustrates how the correlation between the in situ and APP-derived 
albedo decreases with elevation. The reason for the decrease in correlation is because at 
the lower elevation sites the smaller-scale albedo fluctuations are masked by the large 
changes in albedo associated with changes in the snowpack conditions. Thus, there is a 
high correlation between the albedos at these stations. At the higher elevation sites where 
little (e.g. Crawford Point) or no (e g. Summit) melting occurs and the surface albedo 
variations remain small, the fluctuations in the APP albedo that are not related to changes 
in snowpack conditions become more significant and the agreement between the two 
measurements decreases. 

Figure 7 compares all available coincident clear-sky in situ and APP-derived 
surface albedo at the stations listed in Table 2, excluding stations 4, 10, 1 1 and 15. 

Results show that there is a general negative bias in the APP-derived surface albedos on 
the order of 10% as compared with the in situ measurements. However, if the 4% bias 
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observed in the LI-COR measurements at the Swiss Camp is applicable over the entire 


ice sheet, which may or may not be the case, the differences between the in situ and 
satellite measurements would be closer to 6%. 

The overall 10% bias is also apparent in the monthly averaged clear-sky surface 
albedos (Table 3). Again, keeping in mind that the LI-COR albedo are likely biased high, 
the general differences may be somewhat less. In general, the clear-sky monthly mean in 
situ albedos are greater than the satellite-derived albedos at all the stations and at all 
times of the year. Exceptions occur during March and April 1997 at the Swiss Camp, the 
reasons for which were previously discussed; April-June 1998 at NASA-U; and March 
1998 at CP2. At NASA-U and CP2 however, these differences are near zero. 

Although the two measurements typically differ by 10%, there are times when the 
differences are nearly zero, such as during 1998 at the NASA-U station. Similarly, 
differences may greatly exceed 10%, such as during July 1998 at the JAR site, where the 
clear-sky monthly mean albedos differ by 38%. More variability is observed in the 
satellite-derived estimates of the clear-sky monthly albedo as noticed by the larger 
standard deviations in the mean albedo estimates. However, as seen in Figure 8, there 
doesn’t appear to be any clear spatial or seasonal pattern in the magnitude of the 
differences in albedos or in the variability of the APP mean albedo. 
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The general seasonal trends in the monthly-averaged clear sky albedos derived 
from AVHRR agree with the ground-based variations in albedo (e g. Figure 8). 

Typically, the albedo is higher in the spring and fall, and lower in the summer as a result 
of an increase in snow grain size associated with an aging snowpack or snowmelt. At 
most sites, precipitation occurs in late summer/early fall, and thus the snow albedo is 
found to be less during summer as the snow pack ages and before accumulation of new 
snow occurs. At some sites however, such as at Humboldt or the Summit, there is very 
little change in the surface albedo during summer. At Humboldt for example, the 
frequency of snowfall and little melting (up-slope advection from the west results in high 
accumulation at GITS) in this region during summer explains why the summer albedos 
remain high and change very little. 

(5) Discussion 

Comparisons between Greenland AWS station data and the APP-derived surface 
albedo show that the satellite values are on average about 10% less that the in situ 
measurements. This discrepancy is significant in terms of the amount of energy available 
for absorption by the surface. For example, an albedo difference of 0.85 to 0.75 can 
increase the amount of absorbed energy from 1 5% to 25%. The greater amount of 
absorbed energy can subsequently impact surface temperatures and ablation. Thus, 
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accurate parameterization of the snow albedo is important when modeling the surface 
climate of ice sheets. 

Reasons for the differences between the two measurements essentially fall into two 
different categories: (1) measurement error and (2) algorithm error. Measurement errors 
affect both the in situ and the satellite measurements. Comparison of the LI-COR 
measured surface albedo with that from the broadband PSP indicates a positive bias of 
nearly 0.041. This bias partially explains the differences between the ground-based and 
satellite estimates of the surface albedo. Similarly, calibration of the AVHRR instrument 
may also lead to errors in the satellite-derived estimates of the surface albedo. 

Calibration of the AVHRR instrument using the revised post-launch calibration 
coefficients of Rao and Chen (1999) yield albedo values within 5% (relative) to those 
presented in Rao and Chen (1996). Absolute calibration for channel 1 reveals a 3% 
(relative) accuracy as compared with those by Rao and Chen (1999) and Rao and Chen 
(1996). At this time, absolute calibration accuracy for channel 2 is unknown. 

Furthermore, differences in footprint sizes between the two different observations also 
contribute to the differences observed between the satellite and in situ measurements. 

The impact of the large satellite footprint size relative to the ground based measurement 
on the total albedo error is expected to be greatest in regions that experience summer 
melting and/or lake formation, particularly during melt onset, and least at higher 
elevations where little change in surface conditions occur. 
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Errors in the APP algorithm used to derive surface albedo from satellite 
measurements also may lead to errors in the final surface albedo estimate. It was 
suggested previously that some albedo error might be caused by the conversion from a 
directionally dependent reflectance to hemispherical albedo. In the initial processing of 
the APP surface albedo product, the conversion factors of Lindsay and Rothrock (1994) 
were used. Figure 5 shows that with this angular model, a strong dependency of the 
albedo with viewing angle remains, particularly for relative azimuths ranging from 0 to 
120°. This led to extremely large fluctuations in the surface albedo with satellite viewing 
angle. With the Suttles et at. (1988) model, the model used in the APP processing, the 
dependence on viewing angle is reduced, although a slight dependence does remain. 

Both the Lindsay and Rothrock (1994) and Suttles et al. (1988) models use the 
ERBE data. Uncertainties in the ERBE angular correction may lead to an absolute error 
of nearly 2% in the global mean albedo, but the actual albedo error could be much larger 
(Ye, 1993). Recently, Li (1996) showed that problems still remain with the ERBE ADM 
over snow/ice-covered surfaces, particularly during summer. Further improvements in 
angular models are expected with the Clouds and the Earth’s Radiant Energy System 
(CERES) instrument. In this paper, we do not attempt to generate a new set of angular 
models. Owing to the small number of samples used here, we can only provide a 
qualitative assessment of how well the angular models used in the Polar Pathfinder 
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processing work over Greenland. However, with the processing of more of the AVHRR 
data, the possibility will exist to generate a new set of angular models. 

Another source of error related to the directional reflectance of snow, is the 
surface roughness. Recently, Warren and Brandt (1998) showed that the angular pattern 
of sunlight reflected by snow is altered by surface roughness. Roughness in the snow 
surface causes variations in the incidence angle of the sunlight and may even cause 
shadowing. These effects were found to greatest at large viewing zenith angles, 
suggesting that remote sensing of albedo can be carried out accurately without knowledge 
of surface roughness heights and orientation if near-nadir views are used. In the APP 
processing, however, derivation of the surface albedo is not limited to near-nadir views, 
and thus surface roughness is expected to introduce error in the albedo estimates. The 
magnitude of this error however is unclear. Further work is needed to assess the impact 
of surface roughness on estimates of snow albedo from satellite. 

Other errors in the APP-derived surface albedo may be caused by the narrow-to- 
broadband conversion. The albedo at the TOA is defined as the fraction of incoming 
solar irradiance reflected back to space from the earth’s surface and the atmosphere. 
Thus, it is intrinsically dependent upon atmospheric conditions. In the conversion from 
the two narrowband satellite measurements to a broadband albedo, a linear relationship is 
used. However, the conversion factors that linearly relate the two narrowband 
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measurements to the broadband albedo are dependent upon atmospheric conditions. This 


dependence is shown in Figure 9. Aerosol and ozone optical depths affect the conversion 
factor for AVHRR channel 1, and atmospheric water vapor affects the conversion factor 
for AVHRR channel 2. Dependence on solar zenith angle is apparent in both channels. 
Since the conversion from narrow-to-broadband albedo does depend on atmospheric 
conditions and atmospheric path length, a broadband albedo derived under a specific 
atmospheric condition would not necessarily be valid under different atmospheric 
conditions. The error however is predicted to be small at low solar zenith angles (- 0.05 
absolute error in albedo for 0 << 65 °) and large at high solar zenith angles (- 0.15 at 
0 = 80 °). 

Finally, errors may also result in the atmospheric correction. In the atmospheric 
correction, all corrections are made assuming the surface is at sea level. As mentioned 
earlier, the atmospheric correction is performed assuming a constant aerosol optical depth 
of 0.06 and constant ozone amount of 0.24 cm at all elevations and at all times of the 
year. Since the atmospheric correction is performed assuming all surfaces at sea level, 
this method would bias the surface albedos high at the higher elevations. Running a 
simple test using the 6S radiative transfer model ( Tame et al., 1992) for a broadband top- 
of-the-atmosphere albedo value of 0.8, and using the same values for column amount of 
atmospheric ozone, water vapor and aerosols, we find that there would be a 6% 
difference in the model-derived surface albedo between sea level (0.97) and the summit 
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(0.91). Most of this difference however is caused by atmospheric water vapor The 
atmospheric water vapor derived from the AVHRR brightness temperatures at 1 1 and 12 
pm results in greater column amounts at lower elevations than at the higher elevations 
(Figure 10). In this regard, elevational affects are somewhat taken into account. The 
computed precipitable water vapor is also found to vaiy seasonally, with higher amounts 
during summer, so that there is a slight seasonal dependence of the atmospheric 
constituents accounted for in the atmospheric correction. 

Finally, the possibility exists that some residual cloud effects might be present. 
Conditions are assumed clear based on the approach by Box (1997). It is possible that 
thin cirrus or diamond dust are not detected with this method. This would affect the 
instantaneous observations, but would have a much smaller effect on the overall mean 
bias reported here. 


(6) Conclusions 

The surface energy balance of the polar ice sheets is mainly dependent upon the 
radiative fluxes, and during daylight periods, it is primarily dependent on the surface 
albedo. The surface albedo controls the amount of solar energy available for absorption 
by the snowpack and is thus closely linked with surface melt. Since small changes in the 
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surface albedo can substantially alter the amount of solar energy absorbed by the surface, 
accurate parameterization of the surface albedo is essential. 

With the completion of the AVHRR Polar Pathfinder Project, we now have gridded 
maps of the surface albedo for the Arctic and Antarctic spanning the past 18 years. Thus, 
data are now available to support studies of high-latitude climate-change processes, such 
as trying to gain a better understanding the complex relationships between melt and 
accumulation, the energy balance and atmospheric flow patterns. This paper examined 
the accuracy of the APP surface albedo product over Greenland. Results show that the 
APP-derived surface albedo tends to be on average about 10% less than the in situ 
measurements, but may be reduced to something on the order of 6% considering that the 
ground-based measurements are biased high. The differences between the two 
measurements do not appear to have any spatial or temporal pattern. 

Comparisons with ground-based measurements reveal better agreement between 
the two measurements at sites that have a large range in albedo, sites where seasonal 
changes in surface conditions occur. At the higher elevation sites where the surface 
albedo remains relatively constant, the measurement uncertainty exceeds the natural 
variability. This is a result of combined uncertainties in the calibration, atmospheric 
conditions, and the various aspects of the retrieval algorithm. The use of an improved 
angular model to convert the directional reflectance measurement to hemispherical 
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albedo would likely improve the consistency of the albedo measurements. Data from 
multiangle observations, such as from the Multiangle Imaging SpectroRadiometer 
(MISR) or from the Polarization and Directionality of the Earth’s Reflectances 
(POLDER) instruments should lead to improved angular models, as well as improved 
atmospheric characterization and hence improved albedo accuracy. 

Even though the absolute accuracy of the APP-derived surface albedo still falls 
short of the accuracy required for climate models, the general temporal and spatial 
changes in surface albedo are reproduced quite well with the satellite estimates. Thus, 
the present results are useful to GCM modelers because the AVHRR-derived albedos 
provide better estimates of the albedo than what is available otherwise. In conclusion, the 
APP-derived surface albedos provide temporally and spatially consistent estimates of the 
surface albedo over Greenland. Further work may be focused on improving the absolute 
accuracy of the satellite-derived surface albedo. 
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Figure Captions 


Figure 1. Map of the Greenland ice sheet showing locations of the AWS stations 

Figure 2. AVHRR Polar Pathfinder Surface Albedo Processing Steps. 

Figure 2(a)-(d). Surface albedo for the afternoon satellite pass during 1997 at 4 stations: 
(a) JAR, (b) Swiss Camp, (c) Crawford Point, and (4) Summit. Data are shown for all 
coincident ground-based and satellite measurements. No cloud filtering was performed 
for this Figure. 

Figure 4. Surface albedo (bold) and satellite elevation angle at the Summit for the 
afternoon satellite pass during 1997. 

Figure 5(a)-(b). Dependence of the APP-derived surface albedo on satellite elevation 
and relative azimuth angle during 1997 for (a) Lindsay and Rothrock (1994) and (b) 
Suttles et al. (1988) angular models. Results are shown for solar zenith angles ranging 
from 60-70° and three different azimuth angle bins: 0-60° (star), 60-120° (cross) and 120- 
180° (triangle). 
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Figure 6(a)-(d). Comparison between APP-derived and in situ measured surface albedo 
for the afternoon satellite pass during 1997 at 4 stations: (a) JAR, (b) Swiss Camp, (c) 
Crawford Point and (d) Summit. Plotting criteria include solar zenith angles less than 75° 
and t>0.8. 

Figure 7. Comparison between APP-derived and in situ measured surface albedo for all 
coincident clear-sky (x>0.8) cases during 1997 and 1998. 


Figure 8(a)-(d). Monthly-averaged APP and in situ surface albedo for the afternoon 
satellite pass during 1997 at 4 stations: (a) JAR, (b) Swiss Camp, (c) Crawford Point and 
(d) Summit. Monthly means are computed for all data where the solar zenith angle is less 
than 75° and where x>0.8. 

Figure 9(a)-(b). The fraction of solar radiation within the AVHRR bands 1 (a) and 2 (b) 
as a function of changes in solar zenith angle and changes in atmospheric conditions. 


Figure 10. Precipitable water (cm) derived from the AVHRR brightness temperatures at 
1 1 and 12 pm for the afternoon satellite pass during 1997 at 4 stations: (a) JAR, (b) Swiss 
Camp, (c) Crawford Point and (d) Summit, 
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Table 2. Correlation coefficients between in situ and APP-derived surface albedo for 
clear sky conditions, defined as having effective cloud opacity values (t) greater than 0.8 
and solar zenith angles less than 75°. Only results with at least 10 clear sky observations 
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Table 3. Mean clear-sky monthly averaged surface albedo and standard deviations for 
the afternoon pass during 1997 and 1998. Monthly means are computed for all “good” 
measurements where the solar zenith angle is less than 75° and where the effective cloud 


opacity (x) is greater than 0.8. Periods of obvious bad in situ data, such as due to 
instrument riming or instrument tilting, are excluded from the averaging. Number of 


observations used in t 

he averaging is given in parenthesis. 
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AVHRR Polar Pathfinder Albedo Processing Steps 



Figure 2 
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